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The behavior of two magnetic impurities coupled to corre- 
lated electrons in one dimension is studied using the DMRG 
technique for several fillings. On-site Coulomb interactions 
among the electrons lead to a small Kondo screening cloud 
and an overall suppression of magnetic order. For arbitrary 
electronic correlations and large inter-impurity distances R, 
we find a 1/R^ decay of magnetic correlations. 

PACS numbers: 72.10.Fk, 72.15.Qm 



The almost complete understanding of the behavior 
of magnetic impurities interacting with conduction elec- 
trons in metals is one of the major achievements of 
modern condensed matter theory. The impurity mag- 
netism is mainly affected|-|by two distinct mechanisms, 
namely Kondo scrceningp and the. Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction.B While the RKKY 
interaction typically favors magnetic ordering of the im- 
purities, the Kondo effect works against order since it 
tends to quench individual impurity spins. This competi- 
tion characterizes the magnetism in many heavy-fermion 
compounds,EJ where one essentially deals with a Kondo 
lattice of rare-earth impurities carrying a local moment 
due to inner /-electrons. As emphasized by Varma,a the 
basic features of the interplay between the Kondo effect 
and the RKKY interaction are already contained in the 
two-impurity Kondo problem which has attracted much 
attention lately.Q 

So far all studies of the two-impurity Kondo prob- 
lem have ignored correlations among the conduction elec- 
trons. While this is certainly a reasonable assumption 
in most conventional metals, it breaks down if devia- 
tions from Fermi liquid, theory become important. For 
instance, Brugger et alB reported experimental evidence 
for heavy-fermion behavior in the strongly correlated ma- 
terial Nd2_2;Cea;Cu04. This behavior has been modeled 
by considering a lattice of 4f Nd ions interacting with 
strongly cormlated conduction electrons in the copper- 
oxide planes.l3 Furthermore, experimental studies of mag- 
netic impurities in one-dimensional (ID) quantum wires 
(which are described by the strongly correlated Luttinger 
liquid statd3) are. coming into reach due to recent fabri- 
cation advances.l3 

In this paper, as a prototypical example for the ef- 
fects of correlations among conduction electrons, we em- 
ploy the. density matrix renormalization group (DMRG) 
methodH to study the case of ID interacting electrons. 
Recent DMRG studies for the Kondo lattice model have 
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FIG. 1. Two spin-i impurities (large open circles) coupled 
to a ID Hubbard chain (filled circles). For this example, the 
distance of the impurities is _R = 3, and A^ = 12. 



investigated the influence of electronic correlations on 
charge or spin gaps .113 Here we present results for the 
simpler and physically more transparent two-impurity 
Kondo model, with main focus on the impurity spin-spin 
correlations. 

In the absence of spin or lattice instabilities, ID in- 
teracting electrons form a Luttinger liquid characterized 
by a dimensionless interaction strength g. Generally, for 
repulsive interactions, one has < g < 1, and g = 1 is 
the noninteracting (Fermi liquid) value. While the Lut- 
tinger liquid has been thoroughly studied for the clean, 
case and in the presence of elastic potential scattering,lD 
the implications of magnetic impurities are only begin- 
ning to emerge. Oil the one hand, the Kondo effect in 
a Luttinger liquidHI leads to the formation of a ground- 
state many-body singlet. The impurity spin is screened 
by the electrons like in the conventional single-channel 
Kondo effect for noninteracting electrons, albeit with a 
larger Kondo temperature. One finds Tk '~^ (poJ)'^^'"^^^^ 
instead of the uncorrelated result Tk ~ exp{—l/poJ), 
where J > is the antiferromagnetic exchange coupling 
and pq the electronic density of states at the Fermi en- 
ergy. On the other hand, the RKKY interaction becomes 
quasi- long-ranged in a Luttinger liquid. For two spin-i 
impurities described by spin operators-Si and S2, one 
finds the effective RKKY Hamiltonian,ll3 



Hcff — /CS1S2 



where /C - J^R^^ cos(2kFR) ■ (1) 



The 2fcir-oscillatory RKKY indirect exchange coupling K, 
is defined by second-order perturbation theory in the ex- 
change coupling J. Higher-order terms (which are also 
responsible for the Kondo effect) cannot be written in 
the simple form (m). As a function of inter-impurity dis- 
tance R, the RKKY coupling decays slower than the con- 
ventional 1/R Fermi liquid result, which is recovered by 
putting 5 = 1. 

The competition between the RKKY interaction and 
the Kondo effect shows up in the magnetic correlations 
between the impurities, (S1S2). From Eq. (|l|), one would 
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FIG. 2. Numerical results for {S1S2) for J — 1 and several 
U at half- filling. The impurity distance is fixed at _R = A^— 1. 
The inset shows scaled data, where the distance between the 
impurities is measured in units of the lengthscale ^ given in 
Fig. 3 The solid curves are guides to the eye only. 
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3 1 - cxp{K./kBT) 
^ ^ ^' 4 3 + exp(/C/A:i3T) 



(2) 



Due to higher-order contributions in the exchange cou- 
pling J, this expression is expected to saturate for tem- 
peratures below the Kondo temperature Tr-. Therefore 
the ratio between the relevant energy scales, /C/Tff, is 
crucial for the ground-state value of (S1S2). The sign of 
(S1S2) depends on the sign of K, in Eq. (nl), and there- 
fore (S1S2) is also 2fci?-oscillatory. For very small J, 
the RKKY interaction dominates such that (S1S2} takes 
its maximum (singlet or triplet) value. For very large 
J, Kondo screening is effective in quenching the impu- 
rity spins and magnetic order disappears, i.e., (S1S2} -^ 
0. In order to describe the crossover from RKKY 
to Kondo-dominated behavior, a non-perturbative ap- 
proach is mandatory. We note that the unstable fixed 
point separating these behaviors (see Ref. ih is absent 
in the correlated case since the electron-hole symmetry 
required for its presence is broken. 

We have studied two spin-^ impurities S; {i — 1,2) 
coupled to the electron spin density s{xi) at the respec- 
tive impurity location by the standard exchange term 
J s{xi)Si. The ID correlated electrons are modelled as a 
Hubbard chain, with nearest-neighbor hopping t — 1 and 
repulsive on-site Coulomb interaction f/ > 0. The low- 
temperature behavior of the Hubbard chain away from 
half-filling is equivalent to a Luttinger liquid, and the 
appropriate values for g as a function of U and the filling 
factor can be found in Ref. |4[ Directly at half-filling, a 
finite U causes a charge gap leading to a Mott insulator 
instead of the metallic Luttinger liquid state. However, 
since the spin sector remains gapless, the magnetic prop- 




FIG. 3. Lengthscale ^ used in Fig. 2 as a function of U . 
We have chosen ^ = 1 for [7 = 0, and then determined ^ by 
matching data onto the scaling curve. The solid curve is a 
guide to the eye only. 



erties of the Hubbard chain still follow the Luttinger liq- 
uid predictions, but with g = 0. Away from half-filling, 
one has 1/2 < 5 < 1. We have primarily studied the 
half-fillecLcase, for which DMRG is known to be very 
accurate.B To demonstrate that similar conclusions ap- 
ply to the metallic {g > 0) Luttinger liquid, we have also 
investigated the quarter-filled case, for which some data 
are shown below. 

For accuracy reasons, only DMRG results using the 
finite-lattice versioncl and open boundary conditions are 
presented. As the appropriate measure of magnetic or- 
der, the ground-state value of (S1S2) has been calculated. 
We included m ~ 200 states of the reduced density ma- 
trix to study chains of even length N ~ 80, with an accu- 
racy better than 1% for (S1S2) even in unfavorable cases 
(small U and quarter-filling). The two impurities are ar- 
ranged symmetrically around the center of the chain. In 
units of the lattice spacing, their distance R takes odd 
values varying between 1 (impurities at the chain cen- 
ter) and N —1 (impurities at opposite chain ends). This 
geometry is shown in Fig. U The bulk behavior is then 
reproduced in the limit of large N and R ^ N/2. 

In Fig. |, DMRG results for (S1S2) are shown for J = 1 
and various U. Since the data were obtained for half- 
filling and R is odd, RKKY correlations are always an- 
tiferromagnetic such that (S1S2) < 0. We notice several 
features in Fig. ||: (a) Increasing the Coulomb interaction 
U leads to smaller values of |(SiS2)|. (b) The magnetic 
correlations (S1S2) fall off with increasing distance due 
to the decreasing RKKY coupling /C. (c) As shown in 
the inset, all data points can be scaled onto a univer- 
sal curve by using the lengthscale ^ shown in Fig. 0. The 
overall behavior of ^ is quite similar to the Kondo screen- 
ing length vp/Tx, which, roughly speaking, determines 
the extent of the electronic screening cloud around the 
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FIG. 4. (S1S2) as a function of J at half-filling for i? = 5 
and various U . Dotted curves are guides to the eye only. 



impurity. 



Since Tk increases with correlations, the 



observed decrease of ^ with increasing U is in accordance 
with the results of Ref. |ri|. In that sense, electronic cor- 
relations imply a smaller screening cloud. 

As depicted in Fig. ^, DMRG results for the J- 
dependence of {S1S2) at the fixed distance R — 5 ex- 
plicitly demonstrate the competition between RKKY 
and Kondo-dominated behavior. Since almost the same 
curves were found for A^ = 30, finite-size effects seem to 
play only a very minor role here. For small J, the (here 
antiferromagnetic) RKKY interaction leads to a singlet, 
(S1S2) -^ —3/4, while the Kondo effect destroys mag- 
netic order for large J such that (S1S2) — > 0. The qual- 
itative influence of Coulomb correlations can be read off 
from Figs. and B. The magnetic correlations |(SiS2)| 
are significantly reduced by switching on and increas- 
ing the Coulomb interaction U. Regarding the impu- 
rity spins, RKKY-related magnetic order is suppressed 
by electronic correlations. The same qualitative finding 
was obtained for quarter-filling and other values of U or 
R under consideration. 

This behavior can be simply rationalized as follows. 
Switching on the Hubbard-f/ implies an increase in the 
Kondo temperature and enhances the importance of 
Kondo screening. Since the correlated RKKY law (|l|) 
does not change appreciably with increasing U for the 
values of R studied in the DMRG simulations, the in- 
crease in Tk is decisive here and leads to the observed 
partial destruction of magnetic order. 

Next we study the distance dependence of (S1S2) in 
the bulk limit R < N/2, for which DMRG results are 
shown in Fig. |^. For the quarter-filled case and the con- 
figurations considered here (odd R), the impurity correla- 
tions (S1S2) exhibit an oscillatory behavior. We plot here 
only the negative values; similar conclusions are reached 
for the positive ones. For small enough J, the RKKY 
interaction locks the impurity spins in a singlet state ir- 
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FIG. 5. Distance dependence of (S1S2) on a dou- 
ble-logarithmic scale for [7 = 4 and J = 0.25 (circles), J = 0.5 
(squares), J — 1 (diamonds), and J = 2 (triangles). Open 
symbols are for quarter-filling, filled symbols for half-filling. 
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FIG. 6. Power-law exponent a for the data in Fig. 5. 



respective of the distance R. Increasing J causes two 
effects: (a) The Kondo effect leads to smaller values of 
1(8182)1. (b) There is an (approximate) power-law de- 
crease in the distance dependence 



(S1S2) - cos{2kFR)/R° 



(3) 



The exponent a is shown in Fig. ra for the data of Fig. 0. 
Clearly, a — for very small J, and then a continuously 
increases with J. For large J, the limiting value a = 2 is 
approached for arbitrary U and filling factor. 

Since the RKKY coupling A^ — > as i? ^ 00, the large- 
distance behavior should effectively be determined by the 
large- J behavior. Therefore, from our DMRG results we 
expect the asymptotics 



(81S2) - cos{2kFR)/R^ 



(4) 



for any value of J. For J ^ oo, Eq. (Q) holds on all 
lengthscales. Furthermore, additional DMRG results not 
shown here reveal that the asymptotic law (|j) is ap- 
proached faster in the case of strong electronic correla- 
tions. 

The large- J value a — 2 can be understood analyti- 
cally. For J — oo, each impurity spin forms a strongly 
coupled singlet with the conduction electron spin at that 
site (see Fig. ^. This singlet cannot be broken up, and 
therefore the Hubbard chain is effectively cut at the im- 
purity sites. The leading 1/ J contribution to (SiS3)-|Can 
then be computed by open boundary bosonization.ll3 To 
proceed, we write 



(S1S2) = (27r/fcf )2 (s(xi)s(x2)) , 



(5) 



since s(xi) is antiparallel to Si in 1/ J accuracy. We then 
compute {s{x)s{y)) under open boundaries at the impu- 
rity locations (we take Xi = and X2 = R), 

3fc2 

(s{x)siy)} = -^ cos[2fcF(a; ~ y) - f{2x) + /(2y)] 



X [P{2x)P{2y)](^+3)/2 



P{x-y) 



[P{x + y) 



1+9 



- 1 



+ {v ^ -y) 1 

with the functions (we consider kpR ^ 1) 

2^ -1/2 



P{X) = 1 



f{x) — arctan 



2kpR /"Kx\ 



sin(7ra;/i?) 
exp(7r/fci?i?) — cos(7ra;/i?) 



Evaluatina-J;his result for x near xi and y near X2 
reproducesliZI the numerically observed behavior (^. Of 
course, the prefactor in Eq. (m depends on the interac- 
tion strength parameter g. 

To conclude, we have employed DMRG simulations to 
study the two-impurity Kondo problem for interacting 
ID electrons. On-site Coulomb interactions were shown 
to partially destroy magnetic ordering between the impu- 
rities. The main reason for this effect is the increase of the 
Kondo temperature for correlated electrons. For the im- 
purity spin-spin correlations, we obtain a cos{2kpR) / R^ 
behavior at large distances. 

We wish to thank P. Fulde, H. Grabert, H. Schoeller, 
and P. Thalmeier for useful discussions. 
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